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STEADY  VIBRATIONS  OF  WING  OF  CIRCULAR  PLAN  FORM* 

By  N.  E.  Kochin 

The  nonvortical  notion  of  an  Ideal  incompressible  fluid  has  been 
solved  (reference  l)  for  the  case  of  untfora  rectilinear  motion  of  a 
wing  of  circular  plan  form.  The  method  developed  in  reference  1 may 
also  be  generalized  to  the  case  of  the  nonsteady  motion  of  such  wing. 
The  problem  of  the  steady  vibrations  of  s circular  wing  ij  solved 
herein.  The  results  will  be  frequently  referred  to  herein.  The  prob- 
lem of  the  steady  vibrations  of  a circular  wing  was  solved  by  another 
method  by  Th.  Schade  (reference  2). 


1.  Fundamental  equations 

The  wing,  the  motion  of  which  is  under  consideration,  is  assumed, 
as  in  reference  1,  to  be  thin  and  slightly  curved;  its  projection  on 
the  xy-plane  has  the  shape  of  a circle  ABCD  of  radius  a with  cen- 
ter at  the  origin  of  coordinates.  The  principal  motion  of  the  wing  is 
assumed  to  be  a rectilinear  translational  motion  with  constant  velocity 
c parallel  to  the  x-axis.  The  coordinate  axes  are  assumed  as  displaced 
with  the  same  velocity.  On  the  principal  motion  of  the  wing  is  super- 
posed its  additional  harmonic  vibration  of  frequency  o>,  where  the  pos- 
sibility of  deformation  of  the  wing  is  not  excluded.  The  equation  of 
the  surface  of  the  wing  may  then  be  represented  in  the  form: 

z(x,y,t)  = C0(x,y)  + ix(x,y)  cos  cot  + £2(x,y)  sin  cot  (l.l) 

where  the  ratios  t^/ a as  well  as  the  derivatives  3^/dx  and  dC^/dy, 
where  k « 0,1, 2,  are  assumed  small  magnitudes. 

The  fluid  is  assumed  incompressible  and  the  motion  is  assumed  non- 
vortical and  occurring  in  the  absence  of  external  forces.  The  velocity 


*"0b  ustanovivshikhsya  kolebaniyakh  kryla  krugovoi  formy  v plane." 
Prikladnaya  Matematika  i Mekhanika,  Vol.  VI,  1942,  pp.  287-316. 
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potential  will  be  denoted  oy  <p(x,y,z,t)  and  steady  vibrations  of  the 
fluid  will  be  assumed;  that  is,  the  velocity  potential  is  represented 
in  the  form: 


<p(x,y,z,t)  » <j>Q(x,y,z)  + <P1(x,y,z)  cos  oft  + <P2(x,y,z)  sin  a) t 

It  is  evident  that  the  functions  <Pq,  9p  and  satisfy  the 

equations  of  Laplace 


& 


k 


+ 


^k 

Sy2  V 


= 0 


(k  - 0,1,2) 


The  velocity  of  the  particles  of  the  fluid  near  the  leading  edge 
of  the  wing  DAB  is  assumed  to  approach  infinity  as  6'1'  , where  6 
is  the  distance  of  the  particle  from  the  leading  edge,  but  the  velocity 
of  the  fluid  particles  near  the  trailing  edge  of  the  wing  BCD  is 
assumed  as  finite.  From  this  edge  a surface  of  discontinuity  passes 
off  on  which  the  function  <p  undergoes  a discontinuity.  As  in  refer- 
ence 1,  the  problem  will  be  linearized.  Since  the  values  of  the  functions 
and  their  derivatives  are  assumed  to  be  small  quantities  of  the  first 

order,  their  squares  and  products  are  rejected.  The  functions 
9k(x,y,z)  are  further  assumed  tc  have  discontinuities  on  the  infinite 

half- strip  S situated  in  the  xy-plane  in  the  direction  of  the  nega- 
tive x-axis  from  the  rear  semi  circumference  BCD  of  the  circle  S to 
infinity.  The  boundary  conditions  on  the  surface  of  the  wing  a va 
replaced  by  the  conditions  on  the  circle  S located  in  the  xy-plane. 
Everywhere  outside  the  half- strip  S and  the  circle  S the  functions 
<Plc(x,y,z)  are  thus  regular  functions. 

The  boundary  conditions  which  these  functions  satisfy  are  now  set 
up.  On  the  surface  of  discontinuity  3,  the  kinematic  condition 
expressing  the  continuity  of  the  normal  component  of  the  velocity  must 
first  of  all  be  satisfied: 


from  which  is  obtained  the  conditions 


(1.2) 


« 


< 


V 


1 
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The  djnamical  conditions  expressing  the  continuity  of  the  pressure 
in  passing  through  the  surface  of  discontinuity  S are  now  stated. 


If  a stationary  system  of  coordinates  is  employed,  con- 

nected with  the  coordinates  xyz  of  the  moving  system  of  coordinates 


by  the  relations 

x ■ - ct^ 

y = 

z = *1 

t = t-L 

then  the  pressure  may  be  determined  by  the  following 

— 

formula: 

_ _ d<P  p 

P = ’ p3t^  2 

(59  1 

2 * &D*  * fefl 

+ F(tx) 

0-3) 

Since 

&<p  dtp  c^cp  d<p  d<p  ^ <p  _ c 

a 'St  “ c "5x  ^ “ 3x  5y^  _ oy  "c 

Is.  _ Off 

3Z^  dz 

(1.4) 

the  following  equation  will  apply  in  the  movable  xyz  system: 

9<p  d< 

p - - p 3t  + pc  5 

E . £ 

ic  2 

(£)2*(i$Mi?! 

£2J 

+ 

c? 

(1.5) 

When  small  quantities  of  the  second  order  are  rejected  and  the 
magnitude  F(t)  is  not  dependent  on  the  coordinates, 


P = 


pc 


&P 

ox 


or,  on  account  of  equation  (1.2), 


p(x,y,z,t)  = pc  ^2 


+ l pc  ^ puxp2  i cos  + lpc  + sin 


For  briefness,  the  following  notation  is  introduced: 

cd/c  = k 


(1.6) 

(1.7) 


The  condition  of  continuity  of  the  pressure  on 
three  equations: 


then  leads  to  the 
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4. 


on  S 


(1.8) 


The  condition  on  the  circle  S 1b  now  written.  Equation  (l.l)  in 
the  stationary  system  of  coordinates  has  the  form: 

zi  “ to(xrcti>yi) + ^i(xi‘cti»yi) cos  “*1 + s2'xi-cti»yi) 8tn  °)ti 

Hence,  for  the  normal  component  of  the  velocity  of  the  fluid  parti- 
c’  »s  adjacent  to  the  surface  of  the  wing, 


'IZi  dt0  S52 

= - c -gj — c gj-  cos  art  - sin  ait  - c gj-  sin  a)t  + (o£2  C0B 


The  notations 
* c ^ = ■ cl 


- Zi(x»y)  ■*  c(-§^  + k^])  = 2,s(x,y) 


yield  the  boundary  condition 

(a 


z=0 


ZQ(x,y)  t Z^(x,y)  cos  cot  + Z2(x,y)  sin  cot 


which  must  be  satisfied  on  both  the  upper  and  the  lower  sides  of  the 
circle  S and  which  breaks  own  into  the  three  conditions: 


^z=0  = Zk(X’y)  °“  3 


(k  = 0,1,2) 


(1.9) 


The  presence  of  conditions  (1.2)  and  (1.9)  permits  consideration  of 
the  functions  <pk(x,y,z)  as  odd  functions  of  Z: 

^(xjy, -z)  = - *k(x,y,z) 


t 


1 


t 


V 


(1.10) 
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If  it  is  assumed,  in  particular,  that  z = 0, 

<Plc(x>y»0)  = 0 


(1.11) 


in  the  entire  xy-plone  with  the  exception  of  the  circle  S Bind  the 
half  strip  S on  which  9^  undergoes  a discontinuity. 

The  conditions  (1.8),  because  of  equation  (1.10),  assume  the  form*. 

= o ^ « 0 ^ = 0 on  S (l,12) 

Finally,  the  absence  of  a disturbance  of  the  fluid  far  ahead  of  the 
wing  leads  to  the  evident  conditions  at  infinity: 


lim  = 11m  = 

:-+-^r  X-+-5* 


(1.13) 


x -►+• 


The  problem  of  determining  the  function  <pQ(x,y,z)  satisfying  all 
obtained  conditions  for  this  function. was  considered  in  reference  1. 


The  following  equality  is  set  up: 


$(x,y,z)  = fjtayjz)  + i4>2(x,y,z) 


so  that 


Also , 


^(xjy)  + iC2(x>y)  = £(x>y) 

zl(x,y)  + iZ2(x,y)  = Z(x,y)  = - 
The  shape  of  the  wing  will  be  determined  by  the  equation 
z(x,y,t)  = $0(x,y)  + Re^(x,y)e*lu)t^ 


(1.14) 


i(x,y,z,t)  = 4>0(x,y,z)  + Rej$(x,y;z)e‘la>*)  0-15) 


(1.16) 


(1.17) 


■The  functions  <p(x,y,z)  will  then  be  a harmonic  function,  regular 
in  the  entire  half- space  z>  0 and  satisfying  the  conditions: 
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(tiL,  ■ «*-»>  « S <1-“> 

+ ik$  *.  0 on  S (1.19) 

OX 


following  from  equations  (1.9)  and  ( 1.12).  In  the  entire  remaining 
part  of  the  plane  xy  the  following  condition  must  be  satisfied: 

$(x,y,0)  « 0 (l. 20) 

Moreover,  the  following  conditions  must  be  satisfied  at  infinity: 


lim 
x -*-+*• 


lim 


5S  = 


o 


(1.21) 


which  are  the  boundary  conditions  pf  the  first  derivatives  of  the  func- 
tion <t>(x,y,z)  near  the  rear  semi  circumference  BCD  of  the  circle  S 
and  the  condition  that  near  the  forcard  semi  circumference  DAB  these 
derivatives  may  become  infinite  to  the  order  of 


2.  Fundamental  formulae 

In  reference  1 an  expression  was  constructed,  which  depended  on  an 
arbitrary  function  fQ(x,y),  which  determined  a harmonic  function 

q>0(x,y,z)  satisfying  all  the  conditions  imposed  in  the  preceding  section 


The  functions  K(x,y,z,£,T})  and  G(x;y,z,r)  for  z>0  sire  given  by 
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which  are  harmonic  functions  of  x, y,z  where 

r » a /(x~~  c)2  + (y  - n)2  + *2 
R » - x2  - y2  ~ z2)2  + 4a2z2 

In  order  to  satisfy  boundary  condition  (1.9) 

c Zo(*#y)  on  s 

it  is  necessary  to  take 

*0(x*y)  a‘  “ zc^X,y)  + *0^)  (2*5? 

where  g^(y)  is  determined  from  a Fredholm  integral  equation  of  the 
second  kind. 

The  solution  of  the  more  general  problem  of  cteady  vibrations  may 
be  presented  in  a similar  form. 

Thus,  f^(*,y)  and  f.r,(x,y)  denote  two  arbitrary  real  functions, 

continuous,  together  with  their  partial  derivatives  of  the  first  and 
second  order,  in  the  entire  circle  S; 

fjfoy)  + J-f2(x'y)  = f(x'y)  (2-s) 

It  will  now  be  shown  that  the  function 


(2.3) 


(2.4) 


satisfies  all  the  conditions  of  the  preceding  section  except 
condition  (1.18). 


8 


NACA  TM  1324 


The  function  G(x,y,z,r ) * as  shown  by  equation  (2.2),  is  harmonic; 
hence  the  function 


L(x,y,z)  «.  e'ilcxf>  e ***  G(x,y,s)  dx 


J+m 


will  be  a harmonic  function.  In  fact, 


d2L  d2L  ^ d2L  c>G 

dx2  dy2  dz2  ^ 


dx 


= - ikG  - 

dx 


dx 


When  this  egression  is  Integrated  by  parts,  it  is  easily  shown 
that  AL  = 0,  since  both  G and  do  /dx  approach  zero  for  x •. 


It  then  follows  that  the  function  $(x,y,z)  likewise  satisfies  the 
Laplace  equation 


A4>  = 0 (2,8) 

where  from  the  form  of  equation  (2.7)  it  is  seen  that  4>(x,y,z)  is 
regular  everywhere  outside  the  circle  S and  half  strip  3.  In 
exactly  the  same  way  it  is  shown  that  the  conditions  at  infinity  (.1,21) 
and  condition  (1.20)  are  satisfied. 

Furthermore , 


||  + ik$ 


2n 


+ ikK  + 


G(x,y,z,r)  A /a2  - K2  - r\2  cos  y dr 
(£2  + tj2  + a2  - L’a£  cos  ? - 2a7}  sin  y) 


V as  dn 


(2.9) 
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It  is  clear  that  if  x2  + y2>a2  then 


+ ik$  » 0 for  z = 0 


eo  that  condition  (1.19)  is  likewise  satisfied. 

It  thus  remains  to  check  the  finiteness  of  the  derivatives  of  the 
function  $(x,y,z)  at  the  points  of  the  semicircumference  BCD  of  the 
circle  S and  to  establish  the  behavior  of  these  derivatives  near  the 
forward  semicircumference  DAB.  But  near  the  forward  semicircumference, 
the  inside  integral  in  formula  (2.7)  evidently  remains  bounded,  as  do 
its  partial  derivatives;  since  the  first  derivatives  of  the  integral 


J 7 


f(5,r))  K(x,y,z,£,T|)  <15  dr) 


as  established  in  reference  1,  and  a 3 will  again  be  proven,  have  near 
the  contour  of  the  circle  S the  order  5'1/2  (where  6 is  the  dis- 
tance of  a point  to  the  contour  ABCD  of  the  circle  s),  it  is  clear 
that  the  first  derivatives  of  the  function  $(x,y,z)  also  have  the 
order  6~l/2  near  the  forward  semicircumference  DAB  of  the  circle  S. 


For  determining  the  behavior  of  the  function  4>(x,y,z)  ne;-r  the 
rear  semicircumference  BCD,  the  right  side  of  equation  (2.9)  is  trans- 
formed. Denoting  it  by  M(x,y,z)  and  making  use  of  formula  (2.1l)  of 
reference  1 and  the  formula  of  integration  by  parts  (2.14)  of 
reference  1, 


c,y,*)=^  JT  ffc.nMK- 

s 


(42+qS+a2-2a£  cos  y-2an  sin  r)[  2n  J J ^ 


1 

V>~  - n 

Cd 


jL*Ldn  i rrii 
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It  is  evident  that  this  function  remains  finite  near  the  rear 
semi  circumference  BCD.  But  when  the  following  equation  is  integrated, 

+ ikt  ^ M(x,y,z)  (2.12) 

there  is  obtained 


4>(x,y,z) 


eil£X  M(x,y,z)  dx  + 4>(0,y,z)  e*ikx 


(2.13) 


whence  it  is  clear  that  both  the  function  $ and  its  derivative  with 
respect  to  x remain  finite  near  the  rear  semi  circumference  BCD.  The 
derivatives  of  M with  respect  to  y and  z will  be  of  the 
order  6”1/2  near  BCD,  as  follows  from  a consideration  analogous  to 
that  which  was  adduced  previously  for  determining  the  behavior  of  the 
functions  (x,y,z)  near  the  forward  edge  of  the  wing  DAB.  Since 


5? 


eikx  ^ e-ikx 

ay  dy 


it  is  clear  that  the  derivative  d<!>/dy,  and  similarly  <ft/dz,  remain 
finite  near  the  rear  edge  of  the  wing  ECD. 


The  function  (2.7)  thus  satisfies  all  the  imposed  ccnditons.  The 
only  condition  not  utilized  was  condition  (1.18) 


on  S 


(2.14) 


When  the  following  formulas  are  employed: 
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lim 
z -*+0 


") d«  dI> 


- 2«f(x,y) 


lim 
z -*+0 


0 


< aVT 

•^x2  + y2  - a2 


it  is  found  without  difficulty  that  on  S 

(|f)7  o = - f(*>y)  + s(y)  e-ikx 


for  x2  + 
for  x2  + 


where  - 


4\ 


I 


/■ 


' 


«(y)  - 


2^ 


r»2 


-ikx 


(x2 


-fe-2-*  - a2)*1/2  (a2  - ?2  - n2)1/2  co»  r tg.n)  ar  dx  d?  an 


+ yc  + ac 


2 ax  cob  r * 2ay  Bln  r)(C2  + I2  + a2  - 2a?  cos  r - 2arj  Bln 


The  following  equation  is  thus  obtained: 

-f(x,y)  + g(y)  e*lkx  = Z(x,y) 


whence 


f(x,y)  = - Z(x,y)  4-  g(y)  e‘ikx 

Substitution  of  this  value  of  the  function  f(x,y)  in  equa- 
tion (2.16)  yields,  for  the  determination  of  the  function  g(y), 
integral  equation  of  Iredholm 


g(y)  = N(y)  + 1 H(y,q)  g(rj)  dq 


where 


N(y)  = - 


2n3 


glkx  Vs2  - _ t[2  G(x>y,z>r)  cob  r Z(?.n)  dr  dx  d£  di) 

fsjx^  + y2  - afy?2  + TJ2  + a2  - 2a?  cob  y - 2a‘.  sin  r) 


y 2 < p 2 
9 ° 

y > 
(2.15) 


(2.17) 

(2.13) 

an 

(2.19) 

(2.20) 
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with  G(x,y,z,r)  according  to  equations  (2.2)  and 


H(y,n) 


y2  | * 

P ,1k  (x  - C;(x2  + y2  . (t2  £ . n2\^2  eoirardi  dt 

I (it2  ♦ y2  + a2  - 2ax  co*  r - Say  Bin  r)(C^  ♦ 0*  ♦ a2  - 2*1  co*  r - 2*1 

1 * 


■in  r) 


(2.?1) 


3.  Computation  of  forces 

The  pressure  p may  he  determined  from  formula  (1.6),  which  with 
the  notation*  (1-14)  may  he  written  in  the  following  form: 


For  the  computation  of  the  forces  acting  on  the  wing,  it  is  neces- 
sary to  know  the  pressure  on  the  circle  S. 

Because  of  equation  (1.10),  the  pressures  above  and  below  the  wing 
differ  only  in  sign: 


P.  - - P+  (3*2) 

For  clarity,  the  signs  of  the  functions  on  the  wing  will  henceforth 
he  assumed  to  he  the  limiting  values  in  approaching  the  wing  from 
above,  that  is,  for  z«*+  0. 

For  the  lift  force  P the  following  egression  is  obtained 
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(7i 

*0 

a 


But  by  formulas  (2.9}  and  (2.1l),  the  following  equation  applies 
on  the  uoner  side  of  the  circle  S: 


StegSi  ♦ 1 k«(».,y,0) 


K(x»y»o,c,n)  &K  drj  + 


i_ 

2x 


k K(x,y,0,C,rj) 


(3.4) 


V*3  - X2  - y2  Va2  - 42  - q2  cos  r dr 

0?  f y2  + a2  - 2ax  cob  r - 2ay  Bin  r)(5fc  ♦ I2  ♦ a2  - 2a4  cob  r - 2arj  Bin  r) 


d*  drj 


This  ejcpression  Is  Integrated  over  the  entire  area  of  the  cir- 
cle S.  The  order  of  integration  is  interchanged  and  the  two  integrals 
must  be  computed  first  of  all  by  formula  (4.13)  of  reference  1 


A/a2  - x2  - y2 

x2  + y2  + a2  - 2ax  cos  r - 2ay  sin  y 


dx  dy  a 2*a 


(3.5) 


It  will  be  proven  further  that 


SS  - Wa2  - c2  - n2  (3.6) 

S 

For  this  proof,  the  following  function  is  considered: 


K(x,y,z,£,J)) 


d«  dr) 


(3.7) 


Because  of  the  definition  of  the  function  K,  the  function 
F(x,y,z)  is  a harmonic  function  over  the  entire  space  outside  the  cir- 
cle %S.  By  formula  (2.35)  of  reference  1,  the  following  condition  is 
satisfied  on  the  surface  of  this  circle: 

||  = - Zn  on  S (3.8) 


* 
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and  therefore  the  function  (3.7)  is  the  potential  of  the  non vortical 
motion  of  a fluid  corresponding  to  the  translational  motion  of  a cir- 
cular dish  with  velocity  +2 n along  the  negative  z-axis  normal  to  the 
plane  of  the  disk.  This  motion,  however,  belongs  to  those  that  have 
been  studied  in  classical  hydrodynamics,  from  which  can  be  taken  the 
corresponding  expression  of  the  function. 


F(x,y,z) 


K(x,y,z,5,t))  d£  dn 


(3.9) 


\ 


Passing  to  the  limit  z-*  + 0 yields  the  formula 


K(x,y,0,e,t))  d£  dn  = 4 


on  S 


which  is  equivalent  to  equation  (3.6),  since  K(x,y,0,£,q)  is  a symmet- 
rical function  with  respect  to  the  points  M(x,y)  and  N(£,q). 

The * following  formula  is  thus  obtained: 


lit 


a2  - 


coa  r dr _ 

2a  Z cob  x - ain  y 


drj 


S 


C2  + T,2 


(3.10) 
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If  this  egression  and  eimilar  expressions  are  substituted  for  the 
function  <pQ,  obtained  from  equation  (3.10)  for  k = 0,  the  final 

egression  of  the  lift  force  acting  on  the  wing  is  obuined: 


C2  + T]2  + a2  - 


cob  r dr 

2&£  cos  r - 2an  sin  r 


(3.U) 


By  integration  by  parts  and  with  the  aid  of  the  following  formula 


2 n 

P cos  y 

t / £2  + ri2  + a2  - 2a£  cos  y - 


_ 2ng 

2aq  sin  x a(a2  - $2  _ rj2) 


(3.12) 


equation  (3.1l)  may  be  rewritten  in  the  form: 


- VJS’*  2 - - q2|l?e(ikfe"tot)  + 


.!* 


r,  h * F7F7 


cos  y dr 


I" 


a2  - 2a£  cos  r - 2aq  sin  r 


; »d£ 


dTJ 


(3.13) 


In  a similar  manner,  the  formulas  for  the  moments  of  the  forces 
about  the  x-  and  y-axes  are  obtained. 

For  the  moment  of  the  pressure  forces  about  the  x-axis 


■SS- 


M x = / / y(P-  - p+)  dx  dy  = - 2 / / yp+  dx  dy 

S’ 


(3.14) 
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there  is  obtained 


The  order  of  integration  is  interchanged  by  use  of  equation  (3.4). 
It  is  here  necessary  tc  compute  two  integrals.  By  formula  (4.44)  of 
reference  1, 


SS 


— -a^A2--...x-  ZZmm « | naz  Bln  r 

xz  + yz  + a^  - 2ax  cos  x - 2ay  sin  y 5 


(3.16) 


It  will  now  be  shown  that 


y K(x,y,0,e,l)  dx  dy 


(3.17) 


For  this  derivation,  the  following  function  is  considered: ' 


*l(x,y,z) 


K(x,y,z,$.,ri) 


1?  dn 


By  formula  (2.35)  of  reference  1,  the  following  equation  applies 
on  the  circle  S: 


dFi 

ST 


2iry 


(3.16) 


ard  therefore  F1(x,y, z)  Is  the  potential  of  the  motion  of  a fluid  cor- 

responding  to  the  rotation  of  a disk  about  the  x-axis  with  angular 
velocity  -2jt,  a case  studied  in  classical  hydrodynamics: 
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If  7)  K(x,y,z,5,n)  dt) 

s 

2 tjz  y*JR  + a2  - x2  - y2  - s2il  - — — g ^ g — 

j 3(R  + XT  + + ic  + 

a/—  < L*77"7  - E £,*  z -2]  «•»> 


X2) 


Passing  to  the  limit  z-*  + 0 yields  the  formula 

//,  K(x,y,0,C,Tj)  d?  dn  = | yVa2  - x2  - y“  on  S 

s 

equivalent  to  equation  (3.17). 

As  a result,  the  following  formula  is  obtained 


y||  + ik* 


A Jb?  - K?‘  - r\2 


[li+ikf}d« 


dn  - 


*4  ff r i — 5 — fn  ^ s-°°-  r-ar-  - — — 

3n2  II  / ^2  + t]2  r a2  - 2at  cos  r - 2arj  sin  r 

S 1 n 


&z  dr\  (3.20) 


Hence,  for  the  moment  of  the  pressure  forces  about  the  x-axis, 
the  following  expression  is  obtained: 


jjjc  = ~ |££ 


3 n 


//v 


a2  - $2  _ ^ 


^f0 

W 


+ Re 


df 

31 


+ lkf 


a 

2n 


.2  1 

£f0  + Re(e_:UBt  fj 

t 

A sin  r cos  r dr  l 

/ £2  + tj2  + a2  « 2a£  cos  r - 2aq  sin  r 

^ ^ it  — 

dC  drj 
(3.21) 
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or,  on  account  of  the  formula 
2 * 


sin  y cos  r dr 

£2  + q2  + a2  - 2a£  cos  r - 


2n£q 


2an  sin  r a2(a2  - £2  - q2) 


'0 


(3.22) 


the  equivalent  expression 
Mx 


- - f ^SS ^ ' *2  ‘ 1,2  ' n Re(lkfe'li0t)  +f^  [f0  + Re(fe-totj| 


1* 


i. 


sin  r cos  y dy 


S2  + n2  + a2  - 2a£  cos  r - 2aq  sin 


r} 


a e;  (Itj  (3.2?) 


In  the  same  way  for  the  moment  of  the  pressure  forces  about  the 
y-axie 


My  = - ss  x(p_  - P+)  dx  dy  = - 2 SS  xp+  dx  dy  (3. 


24) 


there  is‘  obtained 


My  = 2p=//x[^  + + i k S) 


^dx  dy  (3.25) 


It  is  here  necessary  to  employ  the  formulas 


LS*2  + y2  + 


x*/a2  - x2  - y2  3 „ 4 p 

dx  dy  = - nad  cos  y 

a^  - 2 ax  cos  y - 2ay  sin  y 3 


(3.26) 


(3.27) 


JJ 


x K(x,y,0,£,il)  dx  dy 


8 

3 


2 


As  before , there  is  obtained 


2 * 

ft  cos2  r 

J K2  + tj2  + a‘L  ?a£  c 


2n£‘ 


cos  y - 2aTj  sin  r a2 (a2  - £2  - T]2)  a2 

(3.29) 


also 


%-fr 


Hi 


a2  - e2  - n2  <- 1 f0  - Re 


I*-1*1  (l 


f - ik£  f 


3 * 


a!  rf 
2 n 15° 


+ Re(e-tot  f) 


)1  P -2 2 ^—£°slx 

JJ^  K2  + r\d  + a2  - 2a$  c 


iL 


i cob  y - 2arj  sin  y| 


• iz  dT] 
(3.30) 


The  value  can  now  be  computed  for  the  frontal  resistance  W,  which 
is  composed  of  two  parts.  First,  the  normal  force  (p_  - p+)  dx  dy 

acting  on  an  element  of  th*  wing  dx  dy  will  have  a component  in  the 
direction  of  the  x-axis: 


JIACA  TM  1324 


if 


(p-  - P+)  & dy 


dx  dy 


is  the  equation  of  the  surface  of  the  wing.  Integration  of  this 
expression  gives  the  first  part  of  the  frontal  resistance  in  the  fora: 


In  fact,  the  frontal  resistance  W will  be  less  than  W^,  since 
a suction  force  Wg  appears  because  of  the  presence  of  the  sharp 
leading  edge  of  the  wing  DAB;  therefore. 


W 


(3.32) 


The  suction  force  Vg  is  connected  with  the  presence  of  a strong 
rarefaction  near  the  edge  of  the  wing.  This  rarefaction  is  taken  into 
account  principally  by  the  square  terms  of  the  fundamental  formulas  (1.3) 
or  (l.5)  for  the  pressure  and  it  is  therefore  unnecessary  to  employ 
these  formulas  here. 


The  suction  force  Vg  is  computed  from  the  lav  of  conservation  of 

momentum  applied  to  a thin  filament- like  close  region  t containing 
the  forward  semicircuaference  DAB  cf  the  circle  S;  region  t is 
bounded  outside  by  surface  a and  inside  by  part  S'  of  the  upper  side 
of  circle  S adjacent  to  the  semi. circumference  DAB  and  the  part  S’  * 
of  the  lower  side  of  the  circle  S.  Figure  1 shows  a section  of  these 
surfaces  obtained  by  a passing  plane  through  the  z-nxis. 

Tl^e  equation  expressing  the  momentum  law  is  projected  on  the 
x-exis: 


- V2  - / ip  70 s(n,x)  45 


(J.SS) 
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The  left-hand  side  is  the  sum  of  the  projections  on  the  x-axis  of 
all  the  forces  acting  on  the  volume  of  fluid  considered,  and  on  the 
right-hand  side  is  the  total  derivative  with  respect  to  time  of  the 
component  on  the  x-axis  of  the  momentum  of  this  volume;  this  derivative 
consists  of  two  parts,  a volume  integral  connected  with  the  local  change 
of  velocity  and  a surface  integral  expressing  the  transfer  of  the  momen- 
tum of  the  particles  of  the  fluid  through  the  bounding  surfaces  of  the 
volume  t . 


Equation  (3.33)  may  be  written  both  for  the  stationary  system  of 
coordinates  Oixiyizi  and  for  the  moving  system  of  coordinates  Oxy?.. 

For  the  stationary  system  of  coordinates,  expression  (1.3)  is  used 
for  the  quantity  p;  moreover, 

v*  “If  v“ " Is  <3,M) 


By  the  theorem  cf  Gauss 


(3.35) 


From  equation  (1.3)  and  the  equation  just  derived,  the  following 
expression  is  obtained  from  equation  (3.33)  after  a number  of  simple 
transformations  s 


S'+S"  Sr+S" 


Since  &p/dti  and  dtp/dx  near  the  leading  edge  of  the  wing  are  of 

the  order  6-1/2  ‘ and  bty/Sn  and  cos(n,x)  are  finite  on  the  surface  of 
the  wing,  the  last  integrals  drop  out  when  region  t is  extended  to  the 
line  DAB.  The  following  limiting  equation  is  therefore  applicable: 


For  confutation  of  the  suction  force  W2>  the  expressions  must  be 

found  for  the  components  of  the  velocity  near  the  leading  edge  of  the 
wing  DAB.  The  velocity  of  the  fluid  particles  near  the  leading  edge 
of  the  wing  are  shown  to  he  of  the  order  of  if  5 iB  the  dis- 
tance of  the  particle  to  the  contour  C of  the  circle  S.  From  equa- 
tions (1.15)  and  (2.7)  it  is  evident  that 


9(x,y,z,t) 


r? 


2 71 


K(x,y,z,$,t))  d£  dtj  + 


X(x,y,z,t)  (3.38) 

vhere  the  function  x(x,y,z,t)  and  its  derivatives  remain  finite  near 
the  leading  edge. 

The  behavior  of  the  function  is  now  examined  more  closely 


U(x,y,z) 


■II 


f(On)  K(x,y,z,t,n)  a?  dn 


(3.39) 


near  the  contour  C of  the  circle  S.  Therefore, 

r. 


Su 

s 


H fU,n)  de  ^ 


Since  on  C the  function  K becomes  zero,  the  following  equation 
results 


ff  f(4,n)  ac  an 


-Jh 

UcSJ 


df  , 
SI  d^  d,) 


show! tig  the  finiteness  of  this  integral.  Therefore 


NACA  TM  1324 


s + t{K>^  *k  d>i + °(i) 

S 

where  0(l)  denotes  a magnitude  which  remains  finite  when  5 
approaches  0.  But 


Va^  - x2  - v2  - 


+ R 


dK  dK 


2 


7E 


a 
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The  point  vith  coordinates  (x,y,z)  is  brought  into  correspondence 
vith  the  point  of  the  circumference  C vith  the  coordinates 

Xq  ® a cos  0 y0  a a sin  0 Zq  « 0 


and 


r0*  “ (*()  “ £)2  + (yo  " 0)2  ■ S2  + I2  + a2  “ 2aC  cos  0 “ 2aTl  sin  ® 

(3.43) 

Near  the  contour  C,  the  principal  part  of  the  integral 


Ji(x,y,z)  = — 


2a2  'Va2  - g2  - n2  f(g>n)  dn 


2a2r2  + (a2  - £2  - rj2)(a2  - x2  - y2  - z2  + R) 


(3.44) 


is 


For  this  purpose,  the  following  difference  is  estimated- 

a = Ji(x,y,z)  - N(e) 

The  circle  S is  divided  into  two  parts:  the  circle 

of  radius  a - e;  and  the  ring  Sg  lying  between  the  circumferences  of 
radii  a - e and  a. 
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ro2  - r2  » £2  + q2  + a 2 - 2a£  cos  0 - 2eTj  sin  0 - (x-£)2  - (y-q)2  - z2 

■ 25  cos  a(£  cos  0 -I  q sin  0 - a)  - 62 

2a2(r02  - r2)  - (a2  - £2  - q2)(a2  - x2  - y2  - z2  + R) 

* - 2a6r02  cos  a - (a2  + £2  + q2)  62  - (a2  - £2  - q2)  R 
o 

*jjj  Since 

w 

r02  ^ 4a2  R v 2a5  + 62 

therefore 

|2a2(r02  - r2)  - (a2  - £2  •*  q2)(a2  - x2  - y2  - z2  + R)|<2a6r02  + 2a262  + 

(a2  - <2  - T)2)  R 


Hence  if  |f(£>o)j<M  in  the  circle  S then 


But  by  equation  (2.24)  of  reference  1 


Since 

2a2r2  + (a2  - £2  - q2)(a2  - x2  - y2  - z2  + R; 

= 2a2rQ2  + 2a6r02  cos  a + (a2  + £2  + T)2)  82  + (a2  - £2  - q2)  R 
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hence  for  B < a/ 2 

2a?,r2  + (a2  - £2  - T)2)(a2  - x2  - y2  - z2  + R)  > a2^  (3.46) 


and 


f f a2Va2-g2-T|2  as  dn $ 

J J r02Q2a2r2+(a2l^2-Tj2)(a2-x2-y2-z2+Rj]|  I r0* 

3i  Si 

The  last  integral  evidently  does  not  depend  on  0;  hence  it  may  be 
assumed  that  0=0  and  therefore 


r lxa2 2k  "1 4k&2 2k  + 2n  < 

|3^(a2  - p2)3  . P^jp^Q  3«Y(2a«  - C2)3  ^2a » - *?  3a  3'\^e3 

Similarly 


X/ro2 t*A*  ♦ (?-  <"2  - n2)(i  ^72  - *2  ♦ BO  4 r° 4 ^ ’ **  ’ "2)3  **  4" 


'sir 

and 


2k  a-*  

P A 'tiia?  - p2)3  p dp  d# 

J0  JQ  (a2  - 2ap  cos  d + p2)2  ^ 

a-  e 

/*  2k(s2  + p2)  p dp 

/c  -V(a2  - p2)5 

— p=a-e 

_ 6fla< 

2 J _ **** ...  . . .2 

Jp=o  V2at  • eS  ^ 

vJlcl  ^ 

As  a result,  the  following  inequality  is  obtained 

iAli  5 2nM& + + 
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where  0(a)  denotes  a magnitude,  whose  ratio  to  a remains  finite 
when  8 approaches  zero. 

An  estimate  of  the  second  integral  entering  equation  (3.40)  is 
given:. 


2439 
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hence 


and  for  e = 6 


Uzl 


< 


_2*M  U-  + 

a Va  V\/c  §'  J 


l*T2l 


4nM 

a^a& 


From  equation  (3.40)  and  equation  (3.42),  the  following 
on  account  of  the  estimates  (3.47)  and  (3.49): 

dU  _ K(0)  x - + ocx) 

Jt&R 


In  exactly  the  same  way,  there  is  obtained 


IS  = . H *(«)  y^a2  - It2  - y2  - z2  4~R  , } 

dy  naR  v ' 


Finally, 


But 


« dtj 


dK  __  _ 2z 
3z  jrr3 


arc  tan  A + 


2A 

jt(1  + A2) 


_z_  + z(a2  + x2  + y2  4 
r3  rR(a2  - x2  - y2  - 


where 


ar  V2 


2 


z 


+ R 


Assuming  z>0, 


d£  drj  £ 2* 


(3.49) 
is  obtained 

(3.50) 


(3.51) 
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hence 


J f(M m ten  A+  I ITa2  4 f (C,T|)  d5  dn 


< 2 («  + l)  M 


and  therefore 


du 

35 


*f*f  ■»>  f f . . — . + o(D 

nR^*2  - *z-yz-*2  + R '-'..J  Se2^  + (a2  - £2  - n2)(*2  - x2  - y - *2  + R) 


Again  use  Is  made  of  equations  (3.47)  and  (3.42)  and  the  fact  that 
for  z > 0 


*^/a2  - x2  - - z2  + R 

without  difficulty: 


i a“  + x1"  + y“  + z‘ 


j2  « v2  j.  ,r2  . -2 


jg  _ — t X2  + y_2_t  ^ - R N(e)^/r7J  + x2  + y2  + z2  + 0(l) 

dZ  xRa^Vz" 


(3.52) 


From  what  has  been  said  previously  about  equation  (3.38)  it  is 
evident  that  if 

F(s,t|>t)  ■ fo(S>Tl)  + flU>n)  cos  a»t  + ^z(K>n)  6lourt  (3.53) 


N(9,t) 


F(€j*l>t)'Va2  - g2  - q2  d£  dq 

£2  + ^2  + a2  . 2aC  cos  9 - 2aT)  sin  0 


(3.54) 


the  following  results 
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d®  _ xVa2  - x2  - y2  - z2  + 


SE 


Vi"  n2aR 


yVa2  - x2  - y2  - z2  + 


3y 

d<P  = l 
35  " V"i  it2a2R 


H(e,t)  + o(i) 


N(e,t)  + o(i) 


V2  >t2aR 
(a2  + x2  + y2  + z2  ••  R)  * 


(3.55) 


VR  - a2  + x2  + y2  + z2  N(0,t)  + 0(l) 
or,  in  the  coordinates  6,  0,  a 

^ N(0,t)  cob  0 Bin^|  a) 

S n2^ 

N(0,t)  Bin  0 sln(|  c^ 

^ *2V^6 

d<p  N(0,t)  cos(|  a) 


+ 0(1) 


+ 0(1) 


35 


nz'\]Ze& 


+ 0(1) 


(3.56) 


The  computation  of  the  suction  force  by  equation  (3.37)  is 

considered.  An  arc  D*AB'  of  the  circumference  C is  taken  symmetrical 
with  respect  to  the  x-axis  with  subtending  angle  20q  < n.  For  the  sur- 
face a,  the  part  Oq  is  taken  of  the  surface  determined  by  equa- 
tions (3.4l)  for  constant  0q,  where  0 changes  from  -0q  to  + 0Q  and 
a from  -7t  to  +jt  and  two  bases,*  one  of  which,  g^,  corresponds  to 
0 = 0q  and  the  other,  a2>  corresponds  to  0 = -0q,  where  on  these 
bases  5 varies  from  0 to  6q  and  a from  -n  to  +n. 

On  the  toroidal  surface: 

cos(n,x)  ss  cos  a cos  0 cos(n,y)  = cos  a sin  0 cos(n,z)  = sin  a 

to  to  , s to  to  , x sln(l a) 

' = Sc  C0E(n>x)  + ^ coe(r.,y)  + ^ cos(n,z)  „ — + 0(l) 


-•^2a6 
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Hence  simple  computation  shows  that 


1 

2 


cos(n,x)  dS 


-i-  cos  a cos  6 dfl  da  ♦ 

2s* 


-JL  cos  e sin2  i a W da  + 

2 s*  2 


n2(e,t)  cos  e ds  + o(*\/3jJ) 


In  the  same  manner,  the  integrals  taken  over  the  bases  a±  aqd  03 
have  the  order  0(&q).  Hence  if  6q  approaches  zero,  for  the  suction 

force  developed  along  the  arc  D,AB*,  the  following  expression  is 
obtained 


N2(e,t)  cos  B d Q 


Now  when  0q 
suction  force  Wg 


approaches  n/2,  the  required  expression  for  the 
is  obtained  in  the  following  form: 


w2 


_e_ 

2n3 


n 

N2(e,t)  cos  6 d0 


(3.57) 


The  mean  value  of  the  frontal  resistance  is  found.  Equation  (3.3l) 
shows  that  for  the  mean  value  of  W]_ 
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In  the  same  way,  for  the  mean  value  of  the  suction  force 


1 


where 

(k.  - 0,1,2) 

(3.60) 

(3.61) 

4.  Example 

If  a plane  wing  varies  its  angle  of  attack  periodically  according 
to  the  harmonic  law  so  that  the  equation  of  its  surface  is 

z = (0q  + 0-l  cos  o)t)  x (4.l) 

in  the  notation  of  section  1,  the  following  is  obtained 

Co(x>y)  ■ Po*  £i(*,y)  - 0]*  ^(x,y)  - o 

and  therefore 

Z0(x,y)  «*  - cp0  Zx(x,y)  = - cP-l  Z2(x,y)  = - ckpjx 

(4.2) 

Z(x,y)  = 7a  + 1Z2  = - cp1(l  + i x) 

The  function  f(x,y)  corresponding  to  this  value  of  the  function 
Z(x,y)  is  determined  by  equation  (2.18)  where  g(y)  is  the  solution 
of  Integral  equation  (2.19). 

Consideration  is  restricted  to  the  solution  of  the  inverse  problem 
by  assuming  that 

f0(*,y)  = Ao  f(*,y)  = a + bx 
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where  A and  B are  constant  complex  numbers  and  A$  is  a constant 

real  number  and  the  shape  of  the  wing  is  determined  corresponding  to 
this  function.  By  such  a method  it  i3  possible  to  obtain  aleo  an 
approximate  solution  of  the  direct  problem  of  the  nonsteady  motion  of 
a wing  according  to  the  lav  (4.1)  for  the  case  of  small  frequencies  of 
vibration. 

The  forces  acting  on  the  wing  are  determined.  For  determination  of 
the  lift  force  P,  use  is  made  of  equation  (3.13).  The  following  rela- 
tions are  used 


f/V?  *1  = ! ~3  r /V\P^?-'n2  ^ d„  = 0 

(4 .3) 

as  are  equations  (3.5)  and  (3.26),  yielding  without  difficulty 


P « . lfi£ 
n 


jlte  0 inka3  Ae’^  + 


3 ^ 3 

.2  * 2 71  1 

2 Jaq  + Re(Ae“^Dtj J cos  r dr  + | a3  RejBe-^)^  cos2  r dr| 


or 


p , Aq  + Re(Be-ia>t) 


(4.4) 


The  moment  of  the  pressure  forces  about  the  x-axis  equals  zero  on 
account  of  symmetry; 


MJC=  0 


(4.5) 


If  the  moment  of  the  pressure  forces  about  the  y-axis  is  determined 
by  equation  (3.28)  and,  in  addition  to  tne  previously  mentioned  formulas, 
use  is  made  also  of  the  formula 


SJ 


*a5 
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Oi 

K5 

CM 


*V  *a5  - 


.5* 


I* 


[Ao  + Re(Ae*1<ut)j ^ cos2  r dr  - | a4  Re(Be_laJfc)  L coh^  i r 4f| 


or 

■V  - ^ *0  - ^ «•<**"■*>  - *&£  «*|»"« 

(4.6) 

The  frontal  resistance  is  computed.  First  the  suction  force  is 
computed: 

If 


A = A-^  + iAg  B « + iBg 

according  to  equation  (3.53) 

= Aq  + (A-^  4 B-jO  cos  uyt  4 (A2  + B20  sin  cut 

If  equation  (3.54)  is  applied  and  use  is  made  of  equations  (3.5) 
end  (3.26) , 

4 

N(0,t)  = 2na(A04A1  cos  a)t4  A2  sin  art)  4 j it a2  cos  0 (B^  cos  'i*t+B2  sin  u)t) 


Equation  (3.57)  yields  without  difficulty  the  expression  for  the 
suction  force: 

V2  = 8jt2  b2(Aq  4 A^  cos  0)t  4 A2  sin  cot)2  4 

9 

*2  Jt$  a3(AQ  4 A^  cos  o>t  4 A2  sin  cst)(B^  cos  o)t  4 B2  sin  ost)  4 


■||  n2  a4(B^  cos  cut  4 B2  sin  0)t)2^, 
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or 


w2  - *^|ao2  ♦ \ Ai2  ♦ \ A22  ♦ | aAlBi  + | bA2B2  ♦ ^ *2Bl2  ♦ ^ a2B22  J 
^AoAi  ♦ | «AoB^  cos  ©t  ♦ |2AoA2  + | aAoB  ) sin  art  + 

(l  *12  - ^ *22  + | «*A  - | -A  * If  ‘2®12  - 17  *V) 

(*1*2  ♦ | **lB2  + | "*2*1  + A ‘^’s) 


cob  Burt  + 


sin  Burt; 


(4.7) 


The  total  frontal  resistance  is  obtained  by  the  equation 

V * Wi  - 

where  is  determined  by  equation  (3.3l) 


SJ{ {&  * “[(£  * *“)  * 


= - 2pc 


For  the  mean  value  of  the  frontal  resistance  the  following  is 
obtained: 


W « W-l  - W2 


Refig  e~l  f**  ** 

(4.e) 

(4.9) 


where 


«2  - ^r'  (jo2  + \ Ai2  + | A22  + | aAlBl  + f aA2B2  + ± a2Bl2  ♦ ^ a V) 


WX  = - 2po  /J^^0  + 1 Re| 

s 

For  determination  of  the  functions  6111(1  £(x,y)  character- 

izing the  shape  of  the  wing,  equation  (l .16)  is  used. 
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vhere  by  equation  (2.17)  ia  thie  case 

Zo(x,y)  - Aq  + gQ(y)  Z(x,y)  - - A - Bx  + g(y)e-U«  (4.13) 

and  the  functions  go(y)  and  g(y)  in  this  case  according  to  equa- 
tion (2.16)  have  the  form: 


^ - f? 


X 


.1  A 

(r2  + y2  - a2)  2 (a2  - X2  • n2)2  cos  r <lr  dx  dg  dr> 

(*2  * y2  + *2  . coo  r - 2ay  Bln  r)(t^  + tj2  + a2  ~ 2aC  cos  r - 2*1  sin  r) 


g(y)  - 


2*3 


x 


.1  1 

eikx(x2  + y2  - a2)  2 (a2  - t2  - n2)^  (A  + BC)  cob  r dr  dx  dg  dn 

(x2  + y2  + b2  . 2ax  cob  r - 2*y  ain  r)(C2  ♦ n2  + *2  - 2a?  cos  X - 2arj  tin  r) 


Equations  (3.5)  and  (3.26)  yield 


A/a^-y2  * 


( coa  r dr  dx 

+ y2  - e2  (x2  + y2  + a2  - 2ax  cos  r ~ 2ay  sin  r) 

2 


ra2-y2  | * 


g(y)-S| 


t P.  elkxfA  + 2/3  «B  cos  r)  cos  r dr  d» 

I I , x2  + y2  - a2  (x2  iTy2  + a2  - 2ax  cob  r - 2av  *iu  y) 


► (4.14) 

J 


Integration  of  equations  (4.12)  yields 
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£o(x>y) " c x ■ c x 4 My) 

n \ L 1B\  1 ■ e"1^  iBx  1 
£(x,y)  = ^ iKc 

where  h()(y)  and  h(y)  are  arbitrary  functions  of  y 


kc  ■ c xe’itaC  + h(y)  e 


(4.15) 

-ikx 


The  function  go(y)  was  obtained  in  reference  1,  where,  however, 
errors  slipped  into  the  confutations.  Setting 


y ■ - a cos  0 Hq (0)  * sin  0 gQ(-  a cos  0)  (0<  Q<  ji) 


(4.16) 


gives  in  place  of  equation  (4.22)  of  reference  1 


Hn(e)  =4-8 in  0 


+ -g  sin  6 (in 


V 


] \ 2 

1 + sin  ^ 0^ 

1 - sin  | e; 


+ -g  sin  0 [in 


1 N2 

1 - cos  | 0> 


1 + COS  -g  0J 


1 1 - sin  -g  0 1 1 - cos  *1  0 

cos  -g  9 In j — + sin  ^ 0 In  1 — 


1 + sin  -g  0 


1 + cos  -g  0 


U.17) 


Hence  setting  hQ(y)  « 0 and  Aq  = ac  in  place  of  equation  (4.23) 
of  reference  1 yields 

so<  w>  £■  yiit  - -h  L 


8 if 


^ + y/ 


2*^ 


|ln 

T V2^  - ^ - yj 

ln  a/2?  - \Gr~  _ 'sflL  - ■.<vrrT  \ 

— + + y 2it2  - y ^\f2&  + ^a"+  y j 


/ 


(4.10) 


In  particular  for  y = 0 and  y = + a/ 2 the  following  values  are 
obtained  in  place  cf  those  given  in  reference  1; 
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£(x>0) 


XaH^/T  +1) 


0.9263  ax 


^(X)+-  f) = 

— 1«2(2  ■«V"3)  — ^5  ln23  + -~=.  ln{2+  V^)  + ^5  In  3 - 0.9145  ax 

2jt“  8 n Jl^y  5 n 1 


In  the  same  way,  the  expansion  given  in  reference  1 of  the  function 
HQ(0)  in  a trigonometric  series  in  the  interval  0^0$  n should  be 

replaced  by  the  following: 


Hq(0)  = sin  0 


that  is, 


0-) 


sln(2k  + 1)  0 (l  1 
k{k  + 1)  \3  + 5 


4k 


1+1) 


(4.19) 


H0(e)  = \ 32!c,1  Blr.(?k  f l)  S 

where 

^ 0.9348  P5  = 0.1312  09  = 0.0504 

p3  = 0.2667  (B7  = O.0796  

In  connection  with  this,  corrections  should  also  be  applied  to 
the  numerical  values,  which  are  given  in  reference  1,  of  the  coef- 
ficients Bn  of  the  trigonometric  series  for  the  circulation  obtained 
by  the  usual  theory 


Bj_  = 2.2125  aca 


65  = -0.0296  oca  Bg  = -0.0067  aca 


B3  =5  -0.0934  aca 


By  = -0.0133  aca 


Hence  for  the  lift  force  in  place  of  equation  (4.29)  of*refer- 
ence  1,  the  following  is  obtained: 

P0  = *pca  = 3.4755  pc2a2a 
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which  exceeds  the  accurate  value  by  36  percent.  For  the  induced  drag, 
in  place  of  equation  (4.30)  of  reference  1,  the  following  is  obtained 

w0  • 1.9350  pc2a2a2 

which  exceeds  the  accurate  value  by  87  percent. 

Corrections  are  made  in  the  third  exanjple  given  in  reference  1. 

The  value  of  the  definite  integral  is: 

dy-  ^ - | 1 n2('\/r+  1) 

Hence  in  equation  (4.52)  of  reference  1 the  coefficient  of 
sin  0 cob  6 is  simplified  and  assumes  the  value  -3 n2/®.  In  equa- 
tion (4.53)  the  coefficient  of  sin  20  was  incorrectly  computed; 
its  correct  value  is 

82  = - nr + t = '°-14555 

In  this  connection,  the  value  of  the  coefficient  B2  should  also 
be  corrected: 

Bg  - -0.7436  aca2 

For  the  induced  drag  and  the  moment  of  the  forces  about  the  x-axis, 
in  piece  of  the  values  of  equation  (4.55)  of  reference  1,  the  following 
is  obtained: 

W = 0.4343  pa2c2a4  = 0.5640  pu2c2a- 

the  first  gives  an  error  of  140  percent;  the  second  of  55  percent. 

The  shape  of  the  wing  obtained 

z(x,y,t)  = ^2  x - i gQ(y)  x + 


depends  on  the  frequency  of  the  vibrations  and  is  deformed  during  the 
vibrations.  The  rigid  wing  is  oi  greater  interest. 

It  is  possible  with  the  aid  cf  the  results  obtained  to  obtain  an 
approximate  solution  of  the  problem  of  the  vibrations  of  a plane  cir- 
cular wing  for  small  frequencies  of  vibration. 
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The  case  is  now  considered  of  a wing  varying  its  angle  of  attack 
periodically  according  to  the  harmonic  law  (4.1),  so  that  equation  (4.2) 
holds- 

If 

fo(x>y)  ■ Ao  = a + bx 

equation  (4.2)  yields 

2()(x>y)  ” - Aq  + g(y)  Z(x,y)  - - A - Bx  + g(y)  e*lkx  (4.2l) 


i2-y2  I « 


o0(y)  = S 


cos  r dr  dx 


/x2  + y2  - a2  (x2  + y2  + a2  — 2ax  cos  y - 2ay  sin  y) 


F + * J 2 * 


f p p 3 
a^-y^  -r  « 

A ^ 


Gi(y)  = ^ 


eikx  C0B  y fly  dx 

/x2+  y2  (x2  + y2  + a2  - 2ax  cos  r - 2ay  sin  r) 


/+•  g n 


o2(y)  = 


/a2-y2  | * 


, + -^2  * 


eikx  cos2  r dr  dx 


rx2  + y2  -a2  (x2  + y2  + a2  - 2ax  cos  r - 2ay  sin  rl 

(4.22) 


go(y)  = A0G0(y)  g(y)  = AGx(y)  + BGo(y)  (4.23) 

In  place  of  Gfc(y),  their  mean  values  are  taken  over  the  area  of 
the  wing: 


(Jj^y)  - v2  dy>: 


fa2  - y2  dyV  « / M^W*2  " y2  dy 


(k  - 0,1,2) 
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If  the  frequency  of  the  vibrations  is  assumed  small,  or  more 
accurately,  the  magnitude  ka  is  assumed  small,  the  expansion 

e-ikx  -1  - iioc  - i . 

2 

may  be  limited  to  the  first  two  terms. 

From  equation  (4.2l),  the  following  approximate  expressions  were 
obtained 


Zo(*,y)  - - Ao  + AqgJ) 


Z(x,y)  «-A-Bx+(l-  ikx)(AG^  + BG2) 

Comparison  with  equation  (4.2)  results  in: 

- C0O  » - Aq  + AqGq 


(4.25) 


- = - A + AG^  + BG2 

- c^ik  - - B - ik(A(£  + BG2) 


whence 


Ao 


c0^(l  + 2ikG2 ) 
1 - G^  + ikGg 


cpjikfl  - 2(j±) 

1 - G^  + ikG2 


(4.26) 


The  following  is  confuted 

Aa  . 

Gq  - ~~2  f %(y)  Va2  ’ y2  dy  = \ P G0(-  a cos  0)sin2  G d0 


But  by  equation  (4.16) 


g q(-  a cos  0)  2. 

sin  Q Gq(-  a cos  6)  - sin  Q — — = ttj  Hq(0) 


hence,  expansion  (4.19)  is  used,  yielding 
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Ho(9)  sin  0 «0  - 4 - i)  | - i - ± 


-0.4053  - 0.094 7 Ao  = 1.105cp0  (4.27) 

Equations  (4.26)  shovr  that  In  confuting  G1  it  ie  sufficient  to 

use  the  Jerms  of  first-order  smallness  relative  to  ka,  while  in  com- 
puting G2  it  ie  sufficient  to  use  the  principal  term  not  depending 

on  k.  For  small  ka  the  following  results 

2 In  G2  - G?o  + 0(ka2)  (4.28) 

£&  J 

where  Gn  and  G20  are  the  mean  values  over  the  area  of  the  circle  S 
of  the  functions 


G1  " Go  + ikGn  + 


0^;2a 


n 


and  therefore 


VF 


,v2  -3-  * 

y 9 


Oii(y) 


*IL 


x cob  r dr  dx 


^x2  + y2  - a2  (x2  + y2  + a2  - 2ax  cos  r - 2ay  sin  r) 


(4.29) 


AVa2-y2  j * 


CJ,0(y)  = 2^  ^ r _^^===_= coe*  I 

J J 1 Vx2  + y2 " ® ” (x2  + y2  + 1 


dr  dx 


a2  - 2ax  cos  r - ?ay  si~  r) 


(4.30) 


In  fact, 


°1  - G0  ’ ikGu  = 


G*(x,y >r)  dr  dx  Gy 

(4.31) 


where 

G*(x,y,r)  = 


(eikx  - 1 - ilex)  cos  rVa2  - y£ 


/\fx.2  + y2  - &c  (x2  + y2  + a2  - 2ax  cos  y - 2ay  sin  r) 
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The  Interval  of  Integration  with  respect  to  x is  divided  into  two 
parts:  from^  a2  - to  2a  and  from  2a  to  ®.  Since  for  cl>0 

|eia  - 1 - la|  < a2 

in  the  Interval  Va2  - y2<  x<  2a,  |eikx  - 1 - ikxji  (2ka)2  and 
therefore 


G*(x,y ,r)  dr  dx  dy 


<(2ka)2  Gq<  0.381c2a2 


(4.32) 
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Combining  inequalities  (4.32)  and  (4.33)  yields,  on  account  of 
equation  (4.31), 

| Gi  - C0  - < 0.63a2k2  + 0.12a2k2  In  ~ 

from  which  the  first  of  the  estimates  (4.30)  follows. 

In  an  entirely  analogous  manner,  since,  for  a>0 

| eia  - 1 1<  a 

from  the  inequality 

a ^aS-y2  | x _ 

5,-020  - ^ / P P <«***  - 1)  c°«g  » dx  t 

3»r  / / / ^x2  + y2  - a2  (x2  + y2  + a2  - 2ax  cob  r - 2ay  eln  r) 

2 


the  inequality  is  obtained 


g2  “ °20 


2ka  P P * 'S/a2  - y2*  dy  dx  c_  2ka  f /*  dy  dt 

31,2  ♦ »*  - »z  "3'V.e  Jx 


2to2 

3x 


which  proves  the  correctness  of  the  second  estimate  (4.28). 


The  integral  (4.30)  was  considered  in  reference  1.  The  function 
H-l(0)  of ’reference  1 is  obtained  if 

-a  cos  8)  = H-^e)  (0^0^*) 

For  this  function  the  expression  was  obtained  (equation  (4.36) 
of  reference  1 with  the  correction  of  the  error  appearing  therein) 


3*2 

2a 


sin  0 G 


20 
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The  expansion  of  this  function  in  the  interval  O<0<n  in  a trigo- 
nometric series  has  the  form 


Hl(6)  “ L r2k.+l 

k»0 


sin(2k  +l)0 


where 


rl 


JT 

In  tan  | dx  » -0.69314 


Hence 


^20 


n 

\ I °20^"a  cos  Bln2  0 d0 

^0 


-0.0468a 


(4.35) 


The  mean  value  G^  is  computed.  Integrating  (4.29)  with  respect 
to  r yields 


Gn(y) 


{■ 


_ 

sfx?  + y2  - & ^ 2(x2  + y2)  2(x2  4 y2)  x2  + (y  + a)2 


In 


x(a2 


j£L 


(x2  + y2)(*2  + y2  . tt2) 


-arc  tan- 


/2  - 


2 ax 


4 


dx 


(4.36) 


If 

2 

x = at  y = -a  cos  0 • sin  0 Gn(-a  cos  0)  = H(0) 


then 


ifflBfcJHiiitw  . Jti  ^hggg mm  ~ .fr«irtaUft»  wMriraiiri 
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4/ 


o> 

K5 

-* 

00 


ein  e ( 2 

r 'o)  - / siB  tf  I *t2  _ t COB  0 . t * 4 

| ^t2  - •in2  0 Y 2(t?  + cos2  0)  2(t2  + cos2  0)  t2  + 4 

t2(t2  4-  1 + COB2  0)  t2  - sin2  A 

■7-3 1 — n — rr '4 i — rare  tan >at 

(t2  ♦ cos4*  0)(t2  - sin2  0)  6t 


cos4  ~ 
*“*  \ 


Confutation  of  this  integral  results  in 


H(0) 


- *^sin  " 


sin  - cos 


s\  1 (l  ♦ cos  | sjfl  + sin  | 0) 

!)  +\ ■*  0 InX f4f- K ♦ 

/ A ■ C0B  | 0V1  - *ln  1 


, 1 ♦ sin-5  0 

7 sin  20  in  tan  | - In ~— 

4 " 1 ♦ cos  \ 0 


Further, 


'11 


1556a 


Jt  it 

I / Gll("a  cos  e)  sin2  Q d0  = ■£§  / H(6)  sin  6 d6 

4.  * 4 

The  computation  of  the  last  integral  leads  to  the  result 

°n  ji  * - T + 2 J'  -515-^  d"|  = 1-536  ^ = 0. 

Thus  for  small  ka 

Gj  » 0.0947  + 0.1556ika  G2  = -0.0468a 

Substituting  these  values  in  (4.26)  gives 

a _ «o  1 - 0-0956ika  _ nR  ik(0.8106  - 0.311ika) 

pl  0.9053  - 0.202ika  ~ 0.9053  - 0.“202fka 


fl 

- + 

0 

(+.37) 

(4.38) 

(4.39) 

(4.40) 


48 


RACA  TM  1324 


Thus  for  small  frequencies  of  vibration,  to  a first  approximation: 
A0  = 1.1O5c0o  A - (1.105  + 0.1441ka)  C0X 

B = O.895ikc0^  (4.41) 

For  the  periodic  vibrations  with  small  frequency,  in  accordance 
with  the  law  (4.1)  of  a plane  circular  wing,  the  previously  derived  for- 
mulas may  be  used  for  the  force b where  the  values  A$,  A,  and  B have 

the  values  Just  given.  For  the  lift  force,  the  approximate  expression 
is  obtained  from  equation  (4.4) 


pc.2a2  ^2. 


8130o  + 0^(2.813  cos  o)t  -1.766ka  sin  o)t)>  (4.42) 


) 


The  fluctuation  in  the  lift  force  due  to  the  vibrations  of  the 
wing  thus  leads  the  latter  in  phase,  the  maximum  value  of  the  lift 
force  being  greater  than  the  value  which  was  obtained  in  the  computa- 
tion for  the  steady  motion. 

In  the  same  way,  equation  (4.6)  leads  to  the  following  expression 
for  the  moment  of  the  pressure  forces  about  the  y-axis: 


- pc 


2a3  jl. 47300  + 0i(l.473  cos  cot  + 0.867  kc.  sin  cnt^  (4.43) 

The  component  of  the  frontal  resistance  is  determined  in  the 

given  case  by  the  evident  formula 


a P(0q  + 01  COS  00t) 


that  is, 

Wx  = pc2a2 


^2.8130o2  + 1.4O60J2  + 0O01(5.626  cos  cut  -1.766ka  sin  cot)  + 


1.40601^  cos  2o)t  -O.883012  ha  sin  Scot! 


* 


(4.44) 


The  suction  force  is  obtained  from  equation  (4.7),  restricted  to 
the  first  powers  of  ha, 


W2 


pc2a2  ^1.5540o2  + O.777012  + 0o0l(3.1O7  cos  (nt  + 1.888ka  sin  o:>t)  + 
O.777012  cos  2o)t  + 0. 944ka  0^2  sin  2od^  (4.45) 

The  following  expression  is  obtained  for  the  total  frontal 
resistance: 
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W = - W2  = p-2«2  / 1. 2590o2  + 0.630P-L2  + g0p1(2.519  cos  cut  - 

3.653to  Bin  (Ut)  + 0.6302;,2  cos  2ut  -1.827B’2  to  sin  2u4  (4.46) 

) 


For  the  mean  value  of  the  frontal  resistance 


w = pc2a2  ^1.259p02  + 0.630^!^ 


ihe  flapping  wing  is  considered  such  that 


z = (3qX  + cos  0)t 


In  this  case 


Zo(x,y)  = - cp0 


Z(x,y)  = - lkcp1 


(4.47) 


(4.4-3) 


(4.49) 


Comparison  of  these  expressions  with  equation  (4.25)  shows  that 
in  the  case  considered  it  is  necessary  to  take 


A ----2- 
0 1 - Go 


ikcp1(l  + ikG2) 
1 - Grj^  + ik 


that  is, 


. s\ r—  n „ 1 - ika  0.0468 

Ao  - l-105cp0  A = ikcPl  or9053~~0.~202ito 


k cPl  °1  , 

B = s (4.50) 

1 - + ikG2 


k2cp1(0.0947  + 0.156ika) 
O". 9053  - 0.202ika 


(4.51) 


or,  by  restriction  to  small  terms  of  the  second  order  with  respect  to  k, 


A0  - 1.105cp0  A = ikcp1(l.l05  + 0.195ika)  B = 0.105k2cP1 


(4.52) 


For  the  lift  force 


P = pc2a2  ^2.8130O  + 2.813kpj  sin  a)t  + 0.30Ut2apj_  cos  (4.53) 

and  for  the  moment  of  the  pressure  forces  about  the  y-axis 

My  = - pc2a2  |l.473p0  + 1.473kp1  sin  cut  -0.181k2pp1  cos  ;u^  (4.54) 
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The  component  of  the  frontal  re si stance 

813k3oPi  sin  cut  + O.3Olk2a0o£l  coe 

(4. *55) 


Wj_  a P(3q  = pc2a2 


j2.ei3po2  + 2.< 


The  suction  force  will  be,  with  an  accuracy  up  to  terms  of  the 
second  order  with  respect  to  ka: 


w2  = pa2c2  jl.554p02  + 0.777k2t3L2  -O.376po0ik2a  cob  u)t  + 
3.10”- Vi  sin  cut  cob  2uAJ 

For  the  total  frontal  resistance 

W = pa2c.2  jl.259p02a  -0. 777^0^  -0.294^0^  Bin  cut  + 


2«aR_Q_  r*r»a  /iy4-  -l  H 7771r2fl_  2 nn c. 


O.677k^a0Q£1  cos  cut  7 O.777k^0j^  cos  2ay 
Its  mean  value  will  be 


W - pa2c2 


f 


259p02  -0.777k2p^2 


(4.56) 


(4.57) 


(4.58) 


so  that  a decrease  is  obtained  in  the  frontal  resistance  as  compared 
with  the  wing  which  does  not  execute  a flapping  motion. 
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